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RESEARCH MEMORANDUM 


SOME EFFECTS OF ROLL RATE ON THE LONGITUDINAL STABILITY 
CHARACTERISTICS OF A CRUCIFORM MISSILE CONFIGURATION 
AS DETERMINED FROM FLIGHT TEST FOR A 
MACH NUMBER RANGE OF 1.1 TO 1.8 
By Reginald R. Lund strom and Hal T. Baber, Jr. 


SUMMARY 


A model of a cruciform missile configuration having a low-aspect- 
ratio wing equipped with flap-type controls was flight tested in order 
to determine stability and control characteristics while rolling at 
about 5 radians per second. Comparison is made with results from a 
similar model which rolled at a much lower rate. 

Results showed that, if the ratio of roll rate to natural circular 
frequency in pitch is not greater than about 0.3, the motion following 
a step disturbance in pitch essentially remains in a plane in space. 

The slope of normal-force coefficient against angle of attack Cj^ 

was the same as for the slowly rolling model at 0° control deflection but 
Cjj a was much higher for the faster rolling model at about 5° control 

deflection. The slope of pitching-moment coefficient against angle of 
attack Cn^ as determined from the model period of oscillation was the 

same for both models at 0° control deflection but was lower for the 
faster rolling model at about 5° control deflection. Damping data for 
the faster rolling model showed considerably more scatter than for the 
slowly rolling model. 


INTRODUCTION 


The Pilotless Aircraft Research Division of the Langley Laboratory 
has been conducting a series of free-flight tests through the use of 
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rocket-propelled models to investigate the general aerodynamic charac- 
teristics of a low-aspect -ratio cruciform missile configuration having 
trailing-edge flaps. 

Longitudinal stability, control derivatives, and drag at subsonic 
and supersonic speeds as determined from a similar model which experienced 
very low rates of roll are presented in reference 1. 

This paper presents longitudinal stability and control effectiveness 
for a model which is a replica, as regards configurational geometry and 
mass distribution, of the model for which data have been previously 
presented in reference 1. Since it was considered desirable to check 
the techniques of obtaining longitudinal stability derivatives of rolling 
missiles and to determine the general behavior of the missile when step 
inputs of the pitch control are applied, deflected ailerons were employed 
to roll the model approximately 5 radians per second at supersonic speeds. 
Comparisons are made throughout this report with the data obtained from 
the flight of the model of reference 1. 
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SYMBOLS 

normal acceleration, ft /sec 2 
transverse acceleration, ft/sec 2 

exponential damping constant in e - * 3 ^ per second 
wing mean aerodynamic chord, ft 
body diameter, ft 

acceleration due to gravity, ft/sec 2 
dynamic pressure, lb /ft 2 
body cross-sectional area, 

normal-force coefficient, 
lateral-force coefficient. 


sq. ft 

fn _W_ 

g qA 

ft _w 
g qA 
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resultant-force coefficient corrected for trim, 
( C N " C Rtrim)^ + ( C Y - G Y tr im) 2 ' 


where w m and 
G Ytrim 8X6 ^- e ^ erm i ne ^- using the method of reference 1 

p 

moment of inertia about X-axis, slug-ft 
moment of inertia about Y-axis, slug-ft^ 
moment of inertia about Z-axis, slug-ft^ 

Mach number, v/Vq 
period of oscillation, sec 

total area of rear lifting surfaces in one plane including 
body intercept, sq ft 

trai ling-edge-flap area in one plane, sq ft 
velocity of model, ft /sec 


V C 

speed 

of 

sound in air, ft/sec 

¥ 

.model 

weight, lb 

a 

angle 

of 

attack, deg 

C N t rim 

value 

of 

C]\f when model is at pitch trim point 

^trim 

value 

of 

Cy when model is at yaw trim point 


control deflection, deg 

change in roll angle from the model roll attitude at time 
of previous pitch control step input, deg 

rate of roll, radians/sec 

damped natural frequency in pitch, 2rt/P, radians/sec 


4 


NACA EM L55L16 


:> ■> » > .> • :s 




i i > J .5 '■> 


Derivatives : 
da 


C % = 


per degree 


dC m 

Cm^ = —iii per degree 
oa 


MODEL AND APPARATUS 
Model Description 


A sketch of the model arrangement is presented in figure 1 and a 
photograph of the model in figure 2. The fuselage, which had an overall 
fineness ratio of 12. l6, consisted of a 6.40-inch-diameter cylindrical 
section, and a hoattail rear section. The nose section of the model 
was a 2.25-inch-radius spherical segment and a parabolic section which 
provided a smooth transition from the spherical nose to the cylindrical 
section. The stationary forward lifting surfaces and the rear lifting 
surfaces which will be designated in this report as wings were mounted 
on the fuselage in an inline cruciform arrangement. 

The steel wings of clipped delta plan form were flat plates with a 
thickness ratio of approximately 1.3 percent at the wing-body juncture. 
Leading and trailing edges were beveled with the leading edge being 
swept back 76° 23 ' . Wing panels in the horizontal plane were equipped 
with movable horn-balanced trailing-edge flaps as shown in figure 3 • 
Panels in the vertical plane were identical to those in the horizontal 
plane with the exception that the trailing-edge controls were preset 
as ailerons to a differential deflection of 0° 35 ' to cause the model 
to roll at approximately 5 radians /second at supersonic speeds. 

The flap-type controls in the lift plane, which were connected to 
move as a single unit, were programed in a continuous square-wave pat- 
tern by means of a hydraulic system and motor-driven valve. The two 
control positions were 0° and k. 9 °, measured with respect to the wing 
plane. 


Physical characteristics of the model are presented in the following 
table: 
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W, lb 141. 50 

Center of gravity, rear of station 0 . 46.00 

I X , slug-ft 2 0.312 

Iy, Slug-ft 2 22.28 

I^, slug-ft^ 22.28 

a, ft 0.533 

A, sq. ft 0.223 

S w , sq. ft 3.250 

Sf, sq. ft O.267 

c, ft 2.540 


Instrumentation 

The model was equipped with an NACA eight-channel telemeter which 
transmitted a continuous record of normal, transverse, and longitudinal 
acceleration, angle of attack, rate of roll, control deflection, total 
pressure, and static pressure. Angle of attack was measured by a free- 
floating vane mounted on a sting which protruded from the nose of the 
.model. Rolling velocity was measured by a rate gyro. Total pressure 
was obtained by a total-pressure tube extended from the fuselage ahead 
of the wings and in a plane 45° to the two wing planes. A static-pressure 
orifice was located on the cylindrical section of the fuselage ahead of 
the wings . 

Velocity was measured by a CW Doppler velocimeter and agreed closely 
with that obtained through the use of the total pressure. The model's 
position in space was determined by an NACA modified SCR 584 tracking 
radar set. Atmospheric temperature and pressure were measured by a 
radiosonde which was released immediately after the flight. 


TEST TECHNIQUE 


The model, which was launched from a zero-length mobile launcher 
at a 45° elevation angle, was boosted to supersonic velocity by two 
6-inch-diameter solid-propellant rocket motors which together delivered 
approximately 12,000 pounds of thrust for 3 seconds. After model and 
booster separated, the model was disturbed in pitch by a programed square- 
wave deflection of the trailing-edge flaps. Transient responses to the 
step input of the control surface were continuously recorded in the form 
of time histories as the model decelerated through the Mach number range. 
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PRECISION OF DATA 
Corrections 


Velocity data as obtained by the CW Doppler velocimeter were cor- 
rected for flight-path curvature and wind effect at altitude. The 
magnitudes and directions of these winds were determined by tracking 
the radiosonde balloon. 

In order to obtain the angle of attack at the center of gravity, 
the angle of attack measured at the nose was corrected for model pitching 
velocity by the method of reference 2. Angle-of -attack corrections due 
to combined yaw angle and roll rate were investigated and found to be 
negligible. 


Accuracy 

On the basis of the accuracies of the instrumentation and dynamic 
pressure, the maximum possible errors in M, a, 6, and Cjj are listed 

as incremental values. It should be reiterated here that Cjj is based 
bn body cross-sectional area. 


M 

Limit of accuracy of - 

M 

a 

8 

C N 

1.10 

1.80 

+0.01 

±.02 

+ 0.50 

±.50 

±0.10 

±.10 

±0.30 
+ .09 


These errors, dependent upon telemeter and radar precision, are 
essentially systematic in nature. From a consideration of previous 
experience, probable errors are 50 percent less than those just quoted. 
Parameters dependent upon differences in measured quantities or slopes 
such as are more accurately determined than the previously mentioned 

CG 

errors would indicate. 


RESULTS AND DISCUSSION 

The Reynolds number per foot for this test varied from 12.16 X 10^ 
at M - 1.75* to 6.29 X 10^ at M = 1.10. The atmospheric data as well 
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as the mass and moment of inertia are almost identical for the present 
model and model 2 of reference 1. No data are presented below a Mach 
number of 1.1 because the instrument limits were exceeded, probably due 
to the fact that the model was near roll resonance (0/o> = l) in this 
range. 


The control motion from 5=0° to 5 = 4.9° was very rapid, being 
less than 0.02 second. The roll displacement of the model over this 
time was usually less than 7°. 

Even though the roll rate averages about 5 radians per second the 
wing-tip helix angle is very small. For a roll rate of 5 radians per 
second the wing tip helix angle is about 0.2° at M = 1.1 and 0.1° at 
M = 1.75. 


Time Histories 

The time histories of 6, (%, Cy, and $ as obtained from the 

flight test are shown in figure 4. As may be seen in figure 4, the 
Cn trace is so irregular that no period or damping constant could pos- 
sibly be obtained. Even though the model was disturbed only in pitch, 
quite large and irregular values of side force were induced. Figure 4 
also indicates that the steady-state roll rate at 5=0 is much greater 
than when 8 =4.9°. This was also the case with the much slower rolling 
model of reference 1. This means that either the roll damping was 
increased, because the roll caused a change in the downwash pattern over 
the wing, or the aileron effectiveness was reduced because of the higher 
angle of attack associated with the 4.9° control deflection. In figure 4 
the oscillation on the roll-rate trace immediately following control 
motion from 0° to 4.9° is evidence of rolling moment due to combined 
angle of attack and sideslip and is more noticeable at the lower Mach 
numbers . 

Reduction of the data was carried out by use of the method of refer- 
ence 1 which consisted of plotting Cjj against Cy for each of the 

control pulses and, after accounting for the trim as well as possible, 
developing time histories of C^. Because reference 1 showed that at 

very low roll rates the damped harmonic motion effectively took place 
in the plane in space in which the step disturbance was created, it was 
decided to determine whether that condition applied also to this model which 
rolled at a much higher rate. For this analysis an axis system was used 
similar to that referred to in reference 3 as "pseudo- stability axes" 
for missiles having 90° rotational symmetry. This is the same as a body- 
axis system except that the Y- and Z-axes do not roll with the model. 

For the analysis considered here the position of the XZ plane coincides 
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with the XZ plane of a body-axis system at the instant the control sur- 
face is given the step input but holds this position in space (except 
for translations of the model in which case the whole axis system trans- 
lates) during the analysis of that entire pulse. The roll position of 
the model at any time during a pulse as measured from its roll attitude 
at the beginning of that pulse was determined by integration of the meas- 
ured roll rate. The values of Cr and Cy from the body-axis system 

were then converted to Cr and Cy of the pseudo-stability axis system, 
designated Cm and Cy , respectively, by the following relationships: 

C % = C N cos 0 - Cy sin 0 

Cy^ = Cjj sin 0 + Cy cos 0 


If the damped harmonic motion resulting from the step control input 
remained in a plane, it must be in the XZ plane of the pseudo-stability 
axis system. For this condition if the model had perfect symmetry such 
that it would trim out at Cr = 0, Cy = 0, Cyg would always be zero. 

If, however, the model had asymmetry (such as misalinement or deflected 
control surfaces) a plot of Cyg against 0 would be a sine wave having 


its maxima when the asymmetry is in the XY plane and being zero when the 
asymmetry is in the XZ plane. Sample plots of Cyg against 0 are pre- 
sented in figure 5* The values of and- Cy^. m are the values 

of Cyg at 0 = 90 ° and 0=0°, respectively. The irregularity in 
the curve of Cy^ against 0 is evidence that the actual motion is not 

entirely in the XZ plane. This irregularity became more pronounced as 
the Mach number decreased and thus increased the inaccuracies in deter- 
mining C Ntr . m and Cy tr . m . 


The term Cjjg is made up partly of Cjj due to asymmetry and partly 
of Cjj due to the pitching motion. The value of ACr^, which is that 
part of C^g due to the pitching motion, was determined by the following 
relationship: 

= ( C N - C Ntrim) COS ^ “ ( C Y " C Ytrim) sin $ 

Values of Cr and ACr s against time are shown in figure 6 for three 

typical pulses. The solid curve was obtained by use of the method of 
reference 1 and the circled points were obtained by assuming the response 
to the step input was entirely in the XZ plane of the pseudo-stability 
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axis system. The agreement was good for the first four pulses where 
*/o> was less than about 0.3 • Two of these pulses are shown in fig- 
ure 6(a) and figure 6(h). For the remaining pulses, however, of which 
figure 6(c) is typical, the agreement was rather poor. The greater 
amplitude of oscillation and lower damping constant at the lower Mach 
numbers caused much more energy of the pitching oscillation from the 
previous pulse to remain when the pitch controls moved abruptly to the 
new position. As the Mach number decreased, the roll rate increased 
somewhat so that the model rolled through about 270 ° during each pulse. 
Rolling through 90° or 270° per pulse means that the space position of 
the XT plane (pseudo-stability axis system) for one pulse had about the 
same orientation as the XZ plane from the previous pulse. This prevented 
any possibility of the pitching motion remaining in the XZ plane. This 
could have been avoided by decreasing the frequency of the pulses so 
that the oscillation which resulted from one pulse would have decayed 
to a small amplitude before the next pulse started. However, existing 
theory on rolling missiles such as reference 4 does indicate that at 
these higher values of 0/os the pitching motion created by the step 
input would move appreciably out of the XZ plane even though the pitching 
motions from the previous pulse were completely damped. 


Normal Force Due to Angle of Attack 

Typical plots of C-^ against a are presented in figure 7* "Eh 6 

plots show a lower slope for values of a between 0° and -2° than at 
values of a more negative than -3°» Unpublished wind-tunnel data also 
show this nonlinearity as do the data from the model of reference 1. 

Since p was not measured on this model, it was not possible to make 
plots of Gy against p or Cr against resultant angle. Average slopes 
were measured at a » 0° and a ~ -4° and are presented as symbols in 
figure 8. The solid lines in figure 8 are values of Cjj a obtained from 

the model of reference 1. It should be noted that Cp^ is presented in 
reference 1 instead of the plot of Cjf a presented here. Figure 8 shows 
very good agreement for 0^ between the two models at a « 0° and 
6 = 0°, but considerably higher values of Cw for the present rolling 

model than for the slowly rolling model of reference 1 when a « -4° 
and 5 = k.9°. Inspection of the data from the model reported in refer- 
ence 1 showed that was the same value at a « 0° whether 5=0° 

or 6 = 4.8°. However, as may be seen in figure 8 for the faster rolling 
model at a » 0°, Cpf a is much greater when 8 = k.9° than when 8 = 0°. 

Because the roll rates are greater at the pulses when 8=0° where the 
agreement between the two models is very good it is believed that the 
contribution of the Magnus force coefficients to the total Cpj is small 
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and that neglect of these Magnus forces in reducing the data is not the 
reason for the greater value of Cjf a obtained for the faster rolling 

model. The greater appears therefore to be a direct result of the 

combination of asymmetry (due to the control deflection) and roll rate. 


Pitching Moment Due to Angle of Attack 

The plots of Cr against time shown in figure 6 were typical of 
damped harmonic motion. The damped natural frequency variation with 
Mach number is shown in figure 9» The pitching -moment derivative 
was derived from the faired values of u> and is presented as the symbols 
of figure 10. The solid lines on figure 10 are taken from reference 1. 
The agreement is very good for the 8=0° pulses, but the values of 

for the 8 = 4.9° pulses are usually much less from the model with the 
greater roll rate. The values for the first pulse (8 = 4.9°, 

M = 1.73) appear to be approximately on an extrapolation of the curve 
from reference 1. Figure 4 shows that the roll rate is much lower during 
this pulse than during any other pulse. If the same reasoning as 
for that the Magnus terms are probably small because of the good 

agreement at 8=0° is used, the low values of obtained for the 

rolling model must be due to a combination of asymmetry and roll rate. 


Damping 

The exponential damping constant b is presented in figure 11 as 
a function of Mach number. The solid lines are taken from reference 1. 
Direct comparison of damping constant for the two models is possible 
because of the similarity of the mass characteristics of the two models 
and the similar atmospheric conditions experienced during their flights. 
The test points for the faster rolling model are widely scattered. It 
should be noted however that, if the values of b over a pulse had been 
averaged instead of plotted as individual points, the agreement would 
have looked fairly good. The test points indicate that the damping is 
greater at the greater amplitudes but this may not be true. The tech- 
nique used for determining the damping assumes that the model passes 

through Cm, . and Cy, . at the same instant. When this is not the 
“trim •‘■tram 

case it is virtually impossible to determine the actual trim point, 
especially with the limited number of cycles available at the low Mach 
numbers. This uncertainty of the trim is without doubt a major cause 
of the scatter of test points on figure 11. 
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Control Effectiveness 

The effectiveness of the controls in producing lift is presented 
in figure 12. Because of the variation in normal -force trim over a pulse 
during the present test it was necessary to use methods of determining 
Cjig different from those used in reference 1. One method used was to 

take the value at a = 1° (obtained from plots of Cn against a) 

for 6=0° and 5 = 4.9° and plot it against M, a curve being faired for 
each of the S positions. The difference between the two curves is the 
Cjj due to 5 = 4.9°* The second method involved the use of C N values 

just before and after the control surface changed position with precau- 
tions taken to keep within the frequency response limitations of the 
accelerometer. The AC^/As was then obtained from the relationship: 

= ACjj - Cw Ax, where ACjj and Ax refer to the change in Cjj 

or a between the point before and point after the control surface 
motion. The points were chosen so that Ax was vary small in order to 
obtain better accuracy for AC^/AB. Test points for these two methods 
are shown in figure 12 and the agreement between themis good. Agree- 
ment with values of ACjj/AB from reference 1 is fair and would indicate 

no great variation of ACjj/AB due to roll rate. 

Variation of control-surface pitching effectiveness with Mach number 
is presented in figure 13. This was obtained by multiplying values 
of AHji/AB by the distance from the center of gravity to the control 
surface hinge line in body diameters. No comparison is made with refer- 
ence 1 since the values on figure 13 are merely the values on figure 12 
multiplied by a constant. 
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CONCLUSIONS 


Results from the flight test of a cruciform missile configuration 
which rolled about 5 radians per second when compared with the results 
from a similar model which rolled at a much lower rate indicated: 

1. Stability derivatives may be obtained from a symmetrical rolling 

model where $/ao [ that is, — ] is about 0.3 or less but 

\ Pitch frequency j 

damping data will have considerably more scatter than for a nonrolling 
model. 


2. It would facilitate reduction of stability data if the oscil- 
lations from one disturbance are allowed to decay to a small amplitude 
before the next disturbance occurs. 
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3. The oscillations will essentially take place in a plane in space 
if the input is rapid and the value of 0/o> low. 

4. The value of the normal-force derivative Gw for the faster 

’Co 

rolling model was the. same as for the slower rolling .model at 0° control 
deflection hut was much greater for the faster rolling model at 4.9° con- 
trol deflection. 

5. The value of the pitching -moment derivative as obtained 

from the period of the oscillations was the same for the faster rolling 
model as for the slower rolling model at 0° control deflection but was 
less for the faster rolling model at 4.9° control deflection. 

6 . Roll rate causes no great change in the normal force due to 
elevator deflection. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., December 2, 1955* 
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Figure 4.- Time histories 







Figure 4.- Continued 
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Figure 5.- Variation of lateral-force coefficient, as measured along 
pseudo-stability axis, with roll angle. 





Figure 6.- Sample resultant-force-coefficient time histories 
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Figure 6.- Concluded 
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a, deg 


(a) 5 = 0°. 



a, deg 


a, deg 


(b) 6 = 4.9°. 


Figure 7«- Variation of normal-force coefficient with angle of attack 
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Figure 8.- Variation of normal-force -curve slope with Mach number. 
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Figure 9»- Variation of damped natural frequency in pitch with Mach number 





Figure 10.- Variation of static pitching-moment derivative with Mach number. 
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Figure 11.- Variation of exponential damping constant with Mach number 
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Figure 12. - Variation of normal force per unit control deflection with 

Mach number. 
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NACA - Langley Field, Va. 


Figure 13.- Variation of control-surface pitching effectiveness with Mach 

number . 


